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SECTION  1 
INTRODUCTION 


SCALID  is  a particle  moving  code  that  computes  the  time  dependent 
dynamics  of  the  electron  emission  SGEMP  boundary  layer  in  one  Cartesian 
dimension.  The  calculation  is  self-consistent  in  that  it  is  the  correct 
solution  to  the  nonlinear  equations  of  motion  using  the  local  electric  field 
to  accelerate  the  electrons,  but  it  is  assumed  that  no  skin  currents  flow 
on  the  target  surface. 

The  code  will  accept  as  input  an  arbitrary  fluence,  arbitrary  X-ray 
time  history,  and  an  arbitrary  electron  energy  spectrum  which  can  be  read 
in  either  as  a normal  spectrum  or  as  the  spectrum  assuming  a cos9  angular 
distribution. 

As  the  name  implies,  SCALID  computes  with  the  dimensionless  scaled 
equations  of  motion  which  were  discussed  in  Reference  1.  The  advantage  of 
this  is  that  the  dimensionless  variables  are  independent  of  the  separate 
quantities  fluence,  pulse  width,  material  yield,  etc.  so  that  the  numbers 
computed  and  handled  internally  by  the  code  are  always  of  the  same  order  of 
magnitude.  Consequently  the  user  need  not  decide  on  a different  time  step 
and  a different  cell  size  for  each  problem,  thereby  simplifying  the  use  of 
the  code  and  assuring  a uniform  accuracy  and  resolution. 


C.  L.  Longmire,  and  N.J.  Carron , "Scaling  of  the  Time  Dependent  SC.EMP 
Boundary  Layer",  Mission  Research  Corjioration , MRC-R-262,  DNA  3975T, 
April  1976. 
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The  outputs  are  the  dimensionless  electric  field,  potential,  dipole 
moment,  and  charge  and  current  densities  as  a function  of  position,  and  the 
spatial  integral  of  the  current  density,  all  at  selected  times.  In  addition, 
the  physical  units  of  electric  field,  time,  length,  charge  and  current 
densities,  etc.,  are  printed  out  (ir  both  MKS  and  Gaussian  ESU  units)  so 
that  physically  correct  quantities  can  be  obtained  by  multiplying  the 
dimensionless  output  by  these  units. 

The  results  give  the  electron  emission  current  in  the  absence  of 
target  skin  currents,  i.e.,  in  the  absence  of  any  target  SGEMP  response. 

Hence,  if  the  user  inputs  the  correct  electron  emission  spectrum  from  the 
material  surface  into  vacuum,  the  electron  space  currents  calculated  in  the 
code  are  good  for  a conducting  target  only  if  it  can  be  shown  that  the  target 
response  (such  as  skin  current)  does  not  significantly  alter  the  emission 
electron  dynamics. 

Persons  interested  in  obtaining  the  code  SCALID  should  direct  the 
request  to  the  Defense  Nuclear  Agency,  Washington,  DC  20305,  ATTN:  RAEV/SGEMP. 
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SECTION  2 
BASIC  IDEAS 


The  jihysical  principles  leading  to  the  scaling  of  the  boundary 
ayer  equations  to  obtain  scaling  laws  and  dimensionless  equations  were 
discussed  in  Reference  1.  Here  we  only  sketch  the  relevant  ideas. 

_ -> 

When  X-rays  are  normally  incident  with  a flux  F(t)  (cal  cm  “ sec 
on  a flat  plate  with  material  yield  Y (electrons/cal) , electrons  are 
ejected  at  a rate 


r^j(t)  = Y F(tJ  electrons  cm  “ sec 


(1) 


Let  V be  the  normal  velocity  component  of  the  electrons.  Then  tlie  emission 
angular  distribution  tuid  energy  spectrum  determine  dn/dv  (electrons 
(cm/sec)  ^) , the  distribution  of  v.  The  normalized  normal  velocity 
spectrum  is 


g(v) 


dn/dv 

00 

r dn  , 


(cm/sec) 


such  that  J'g(v)dv  = 1.  The  average  normal  velocity  is 

OO 

V = J V g(v)  dv. 

0 


(T) 


(31 
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Let  the  flue nee  be 


and  1 et 


({)  = j"  F(t)  dt  cal/cm“ 
0 


(4) 


= i)ulse  full  width  at  half  maximum  (FWIIM).  (5) 


be  the  FIVHM  in  seconds.  Then  is  a naturally  occurring  characteristic 

time  in  the  problem,  and 


y<p 


cm 


(6) 


is  a naturally  occurring  number  density.  Similarly, 


>'l  = V cm  (7) 

is  a naturally  occurring  unit  of  length. 

Thus  the  source  specifies  a characteristic  velocity  v , time  Tj^ , 
length  Lj  , and  number  density  . These  quantities  can  serve  as  units 
in  terms  of  which  one  c;m  measure  all  dimensioned  dynamical  variables. 


However,  there  is  another  naturally  occurring  unit  of  time,  the 
plasma  period 


T = 
P 


47T  e N 


sec 


(8) 
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:ind  a corresi)onding  distimcc,  the  Debye  length 

SI,,  = V T cm  . (9) 

D p 

T is  the  time  scale  for  the  emission  rate  to  change,  while  T is  the 
h ^ ’ p 

time  scale  for  the  internal  electron  plasma  dynamics  to  change. 


In  any  particle  moving  code  that  calculates  the  boundary  layei 

structure,  the  time  stej)  At  must  be  small  compared  to  both  and 

In  SCALID,  the  unit  of  time,  LIT,  is  chosen  as  the  smaller  of  T and  T,  , 
p h 


IJT  = M1N(T  , T.  ) sec. 

p h 


(10) 


and  the  time  step  is  a small  fraction  of  UT.  Real  time  is  measured  as  a 
multiple  of  UT,  thereby  defining  dimensionless  time,  denoted  here  by  a 
prime,  t': 


Real  time  = 

UT*t  * 

(11) 

Dimensionless  lime 

= t * = Re  a 1 

time/UT 

(Id) 

Real  time  step 

= At  = small 

fraction  of  111 

(15) 

Dimensionless  time 

step  = At ' = 

small  fraction  of  unity 

(14) 

The  time  in  the  code  is  dimensionless  time  t^  , and  the  time  step  is  At* 
(called  DT  in  the  code)  and  is  small  compared  to  1.  i'ypically,  DT = At  ' = 0. 025 
is  used.  Ui'  is  calculated  in  the  code,  but  DT  is  user  suii])lied.  DT 
usual  !>■  need  not  vary  from  problem  to  problem. 

Since  the  problem  has  only  one  characteristic  velocity,  v , once 
a unit  of  time  has  been  chosen  (UT)  the  unit  of  length  is 

UL  = vUT  cm.  (15) 
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The  imit  of  velocity  is.  of  course, 


ll\  = V 


cm/ sec 


( H)) 


If  the  X-ray  juilse  is  not  sinipl)’  a power  of  time,  'l'|^  exists,  and  the  unit  of 
number  density  UN  is  Nj  of  htjuation  (6), 


UN  = N,  cm 


- .■) 


( 1~) 


If  the  pulse  continues  to  rise  as  a power  of  time,  does  not  exist  and 

UN  is  comiHited  as  in  Section  41). 


The  unit  of  emission  rate  is  therefore  determined  to  lie  N^v 
Hie  dimensionless  time  histors’  I (i)  is  defined  by 


■It)  = f- 
h 


(U 


(1«) 


w h e re 


= t/T. 


119) 


It  is  normalized  to  unit)’. 


dl 


0 


l-’D) 


a))d  the  full  width  at  half  maximum  of 
hquation  (1),  c;in  now  be  written 


I ' ( i ) is  un  it)'.  1 he  cm  i ss  i on  rate 


r^jlt)  = N j V l-'l  t ) 


(-M) 


so  that  l’*(l)  is  the  dimensionless  emission  rate. 
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With  g(v)  the  normalized  spectrum  of  hiiuation  (2),  a dimensionless 
spectrum  c:ui  he  defined  as 

g'tv'l  = V g(v)  , (22) 


so  that 


Jg'(v')dv 


1 , 


125) 


where 


V = 


(24) 


is  the  normal  component  of  velocity  measured  as  a multiple  of  v 


In  a particle  moving  code,  each  particle  on  omission  rejiresents 
the  charge  emitted  by  all  electrons  with  velocity,  say,  between  v'  and 
v'  + Av'  and  during  time  between  t'  ;md  t'  + At'.  The  weight  attributed 
to  such  a particle  is  therefore 


AW  = g'(v')  I''(x)  Av"  At' 


(25) 


with 


h(|uation  (25)  is  the  weight  of  an  emitted  particle. 


(2b) 


Once  emitted,  the  jiarticles  are  moved  by  the  dimensionless 
e<|uations  of  motion  derived  in  Reference  1: 


dv' 
dt ' 


a 


cn 


n'(x',t  ')dx' 


(27) 


(28) 
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Here  n'  is  the  dimensionless  number  densit)’  (number  densit>7Nj).  liquation 

(J8)  is  the  correct  expression  for  the  dimensionless  electric  field  if  the 

luiit  of  time  Hi' = i',  , as  in  Reference  1,  and 

h 


= 


(29) 


In  SCALin,  the  unit  of  time  is  given  by  liiiuation  (10),  and  instead 
of  f.cjuat  ion  (28),  the  dimensionless  field  is  given  by 


where 


a'  = - aJ"  n'(x',t')dx" 

X ' 

mi  . 

T “ 

I’ 


(30) 


(31) 


Ii(luations  (27)  and  (30)  are  the  dimensionless  equations  of  motion  used  in 
the  code. 
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SECTION  3 

GENERAL  CODE  STRUCTURE 

After  all  input  ciutintities  are  read  in,  the  space  grid  is  set  up. 
The  space  grid  is  used  for  calculating  the  electric  field  as  a function  of 
position.  Since  dniamical  tiuantities  are  expected  to  vary  more  rapidly 
near  the  surface  than  farther  away,  the  space  grid  has  logarithmic  sjiacing, 
with  each  cell  a constant  factor  thicker  than  the  previous  one.  If 
man>’  cells  are  used  (one  may  use  uji  to  500  cells),  C must  not  be  much 
larger  thaji  imity.  If  there  are  NC  cells,  the  ratio  of  the  thickness  of 
the  last  cell,  first  cell,  is 


(32) 


If  NC  = 500,  this  will  be  unwieldly  large  if  C is  greater  than  about 

1.005,  for  which  ~ 12. 

NC  1 

Hie  first  cell  nearest  the  surface  will  be  the  thinnest.  The 
first  few  cells  should  not  be  made  too  thin.  Ibey  should  be  small  compared 
to  a Debye  length  (that  is,  small  compared  to  unity  in  the  code's  dimension- 
less variables)  but  large  enough  to  contain  on  the  order  of  100  i)articles 
per  coll  for  good  statistics  when  computing  charge  ;uid  current  densities, 
(llie  number  of  particles  per  cell  is  also  affected,  of  course,  by  the  number 
NV  of  particles  emitted  every  time  stej).  ) We  find  good  results  if  the 
first  cell  thickness  is  -0.05  to  0.1. 


n 


riic  in])ut  spectrum  is  read  in  as  an  enerfjy  spectrum  with  enerjiles 
in  keV.  Hie  injiut  energy  spectrum  may  be  either  that  for  a cos0  emission 
angular  distribution  ( 1 CN  = 0)  , or  it  ma>'  bo  the  sjiectrum  for  the  normal 
component  of  energy  (ICN=  1).  If  the  former,  the  spectrum  is  internal  !>• 
converted  to  a normal  spectrum. 

The  spectrum  is  then  converted  to  a normal  velocity  spectrum 
dn/dv  , and  the  input  energies  are  converted  to  normal  velocities.  Tho 
average  normal  velocity  v and  energy  are  then  computed. 

Ibe  spectrum  and  velocities  are  now  normalized,  velocities  (as 
multiples  of  vj  are  the  dimensionless  particle  emission  velocities.  If 
■MV  energies  were  originally  read  in,  there  are  NV  jiarticles  emitted  every 
time  step.  'Hiey  all  have  weight  l/.W. 

Before  entering  the  time  loop,  the  dimensional  units  of  length, 
time,  velocity,  electric  field  strength,  potential,  charge  densit>',  current 
density,  etc.,  arc  computed  and  printed  out. 

Ibe  time  loop  is  then  entered,  with  the  new  particles  being 
emitted  first.  Tl'.eir  weight  is  determined  by  the  time  history  taken  from 
function  .1-1.  Then  all  particles  arc  moved,  using  old  velocities  ;uid  old 

accelerations.  .\ew  fields  arc  then  calculated  at  cell  edges. 

If  the  cycle  is  a print  cycle,  ciuantities  of  interest  that  are  not 
needed  in  the  y>article  moving  looj)  arc  comjnited,  such  as  current  and  charge 
densities  and  potentials,  and  are  printed. 

Hie  code  structure  is  e.xhibitcd  in  the  flow  diagram  in  i'igure  1. 
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SECTION  4 
INPUT 

Here  we  discuss  the  input  card  format  and  the  way  the  space  grid 
and  emission  velocities  are  computed. 

A.  INPUT  CARD  FORMAT 

Four  basic  cards  are  read  in  per  problem,  plus  an  additional  NV 
cards  specifying  the  electron  emission  spectrum  if  it  differs  from  the 
previous  problem  or  if  this  is  the  first  problem. 

The  four  basic  cards  are: 

CARD  1 

An  80  character  alphanumeric  title. 

Tn'LF.(I),  1 = 1,8 
Format  (8A10) 

CARD  2 

Fight  fixed  point  parameters. 

Format  (1415) 

NT  = number  of  time  cycles  this  problem  is  to  run 
NC  = total  number  of  cells  in  space  grid,  <500 
IR  = 0 - Read  electron  emission  spectrum  on  following 
cards 

= 1 - Use  same  spectrum  as  previous  problem  (IR  Ccin 
equal  1 only  if  this  is  not  the  first  problem) 
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ICN  = 0 - the  electron  energy  spectrum  is  for  a cos9 

angular  distribution 

= 1 - the  electron  energy  spectrum  is  the  normal  energy 

spectrum 

NV  = number  of  spectrum  points  read  in,  <100 

IT  = 1,  2,  3,  ...  = choice  of  time  histories,  see 

FUNCTION  FT 

= 1 - X-ray  pulse  rises  as  a power  (POW)  of  time 

= 2 - symmetric  triangle 

= 5 - asymmetric  triangle.  Rises  to  peak  in  0.5*FWHM, 

falls  to  zero  at  2.0*FW11M 

= 4,  5,  6,  ...  - user  must  define  a time  history  in 

FUNCTION  FT 

IP  = print  cycles.  A printout  will  occur  every  IP  cycles. 

LP  = 0 - this  is  last  problem 

= 1 - another  problem  follows 

CARD  3 

Five  floating  point  variables. 

Format  (8F10.3) 

DT  = dimensionless  time  step  (typically  l)T=;  0.025) 

XMAX  = dimensionless  distance  from  surface  to  end  of  grid. 

Essentially  the  number  of  Debye  lengths  in  the  grid. 
(Typically  20  $ XMAX  5 100) 

XI  = Space  grid  logarithmic  spacing  ratio, ^1.0.  If 

XI  = 1.0  all  cells  are  the  same  size.  If  XI  >1.0  , 
each  cell  is  a factor  XI  thicker  than  the  precceding 
one.  If  NC  is  large,  XI  must  not  be  too  large  for 
practical  purposes. 
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I’OW  = power  of  time  of  X-ray  time  history  (used  oiil>-  if 
IT  = 1) 

I’ROBNO  = problem  number 
CARD  4 

Tlirec  floating  point  variables. 

Format  (8F.  10.3) 

VIH  = material  yield  in  electrons/calorie 

FI.IJ  - If  IT  i 1,  l-I.IJ  = X-ray  fluence  in  calories/cm“ 

If  IT  = 1,  FLU  is  tlic  coefficient  determined 
(by  the  user)  from  the  equation 

X-ray  flux  | ■ ^ FLU  ( t ime (sec ) ) ’’^"  . (33) 

\ cm^  sec / 

F'WIIM  - If  rr  / 1,  FWIIM  = X-ray  i)ulse  full  width  at  half 
maximum  in  seconds. 

If  IT  = 1,  ITVIIM  is  not  used. 

lliese  four  cards  are  read  for  each  problem.  In  addition,  if 
IR=0,  NV  more  cards  are  read  in  whicli  siiccify  the  electron  energy  spectrum 
in  format  2ld0.3.  Hic  first  numlier  on  each  card  is  the  electron  energ\’  in 
keV.  Fhe  second  number  is  the  value  of  the  energy  spectrum  at  tliat  energ)’, 
in  arbitrary  units,  i'he  energy  spectrum  is  normalized  within  the  code,  so 
no  si)cc.ial  units  need  l)e  used  on  input.  The  correct  normalization  informa- 
tion is  contained  in  YIF. 
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B.  SI'ACH  CRin 

The  space  grid  is  set  up  first.  The  thickness  of  the  I-th  cell 
is  called  XB(1).  Tlie  first  cell  is  given  by 

XB(1)  = XMAX  (34) 

(XI  ^-1.) 


and  thereafter. 


XB(1)  = XI*XB(I-1)  , I = 2,  3,  NC.  (35) 


An  array  SL(J)  is  defined  as  the  distance  to  the  J-th  cell,  by 


SL(1)  = 0. 

.1-1 

SL(.J)  = ^ XB(I)  , J = 2,  3,  ...  NC+1.  (36) 

1 = 1 


The  space  grid  is  sketched  in  Figure  2. 

Since  the  grid  is  logarithmically  spaced,  the  bins  fall  at  incon- 
venient locations  for  printing  quantities  as  a function  of  .x.  Instead, 
output  quantities  are  computed  every  0.1  length  units  using  an  array 
KIO(I)  defined  as  follows.  Mark  every  tenth  length  luiit  on  the  grid,  as  in 
Figure  3.  Fach  mark  will  fall  in  some  bin.  K1()(I)  is  the  bin  number  in 
which  the  mark  for  x = 0.1*.J  falls.  For  example,  in  Figure  3,  K10(7)  =6, 
and  K10(2)  =3,  meaning  0.7  length  units  falls  in  bin  6,  etc. 

When  computing  output  quantities,  array  klO  is  used  to  locate 
the  bin  corresponding  to  every  0.1  length  imits.  Output  qu;mtities  are 
then  computed  by  averaging  over  that  bin  and  two  adjacent  bins,  as  discussed 
later  in  Section  6. 
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c:. 


t:MissioN  vr;Locrnr;s  ,v\d  wiiicirrs 


As  read  in,  the  emission  spectrum  is  an  energy  siiectrum  with 
energies  in  keV  and  either  a normal  spectrum  (IC.M=  1)  or  the  energy  spectrum 
for  a cosB  angular  distribution  [ICN  = 0).  If  the  latter,  the  spectrum 

■) 

dn/dt;  is  converted  to  a normal  energy  spectrum,  dn/dw,  where  H=(l/2)mv“  , 
and  w=  (l/2)inv“  , whci'c  is  the  normal  veloeit)'  component,  according 

to 


dji 

dw 


(37) 


Hi  is  normal  energy  spectrum  is  then  converted  to  a velocity 
spectrum  b\’ 


Y 


w 


me 


V 


= cjl  - 1/Y- 


dn 

dv 


mv  Y 


3 dn 
dw 


(3H) 


whore  we  .have  dropi)od  the  subscript  x 

2 2-1/'’ 

relativistic  factor  y = (1-v  /c  ) 


on  v.  Here  y is  tlie'  usual 
-) 

, and  Vi  - (.y-1)  1110“. 


i'he  average  normal  velocity 


00 


0 


(39) 


is  then  computed,  and  the  spectrum  is  normalized  to  unit  integral,  as  in 
liquations  3 and  2.  Velocities  arc  now  normalized  to  v.  nicse  velocities 
arc  called  Vi;(I),  ;md  the  spectrum  CV(I),  1=1, NV. 
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l inall}-  the  spectrum  is  divided  into  bins  of  C(|ual  area,  ;md  the 
velocities  \'h(I)  are  redefined  as  the  central  velocities  of  those  bins. 
Hiis  j;ives  to  each  emitted  particle  the  weight  CV'(  1 ) = 1/NV. 


i).  I’llVSlCAl.  IWITS 


SCALin  calculates  with  dimensionless  ijuantities.  To  get  jibysical 
results  these  numbers  must  be  multiplied  h>'  the  api)ro]iriate  dimensional 
unit.  I'hese  units  are  calculated  and  [irinted  out. 


If  tlie  pulse  does  not  rise  as  a |50wer  of  time  (IT  4 1)  then  I'MIM 
e.xists  and  the  unit  of  number  density  is  as  in  liciuation  (b). 


UN' 


V (cm/sec)  • ITVIIM(sec) 


\ cm-’ 


(40) 


and  the  unit  of  time  is 

UT  = MIN(T^.I'WIIM)  sec, 


(41) 


where 


T 

I> 


sec 


(42) 


is  the  plasma  period. 


If  the  pulse  rises  as  a ])ower  of  time,  t^'^''^,  ( IT  = 1 ) then  there 
is  no  TWllM,  and  the  unit  of  time  is  taken  as^ 

1 


UT  = / ^ \ 

\ Tire*' -Yield- FLU  / 


2+ now 


sec 


(4.3) 
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The  unit  of  number  density  is 


UN  = (UT)’’°"  (^) 

V ' cm'  ' 


In  l-.qiiat  ions  (43)  and  ( 44),  the  quantity  Fldl  is  that  from  liquation  (3 


llie  unit  of  velocity  is  always 


ll\'  = V 


;md  the  unit  of  length  is 


IJL  = IJV*UT 


The  unit  of  charge  density  is 


UR  = 4.8xl0'^^’  UN  ^ 
cm' 


or 


URM  = UN 


in  NIKS  units.  Other  units  arc  computed  similarly. 


(4  4) 

-’0. 

(45) 

(4fa) 

(47) 

(48) 
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SECTION  5 

PULSE  TIME  HISTORY 


'niree  time  histories  arc  built  into  SCALU).  The  user  chooses  one 
b>-  the  value  of  IT  on  input.  iiicy  arc  contained  in  ITINCTIO.V  FT. 

IT  = 1 

ilus  is  for  an  X-ray  pulse  that  rises  in  time  like  a 
power  of  t.  Tlie  power  is  read  in  (I’OIV)  on  card  .T. 
iiic  X-ray  flux  is  given  by  liquation  (?>?>). 

IT  = 2 

A symmetric  triangle  time  history.  As  computed  :uid 
used  in  !■  UNCI  ION  FT  its  Tull  width  at  half  ma.ximum 
is  unity,  and  its  integral  is  also  unity.  ihus  it 
rises  to  value  1 at  time  1 and  falls  to  zero  at 
time  2.  ilie  correct  scale  of  time  is  taken  care  of 
when  FUNCi'ION  FT  is  called,  using  F'WilM  read  in  on 
card  4. 

1 T_=^  .3 

■An  asymmetric  triangle.  In  FUNCTION  FT  its  full  width 
at  half  m.iximuin  is  1 as  is  its  integral.  It  rises 
to  1 at  time  0.5,  and  falls  to  zero  at  time  2.0.  Thus 
the  rise  time  is  1/4  of  the  full  pulse  length. 
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I'ho  user  cliooses  one  of  these  time  histories  by  its  stiape,  and 
reads  in  the  real  full  widtii  at  half  maximum  as  MVIIM  on  card  4. 

I'he  coding  in  lUNCTION  I'T  allows  the  user  to  add  additional  time 
histories  (14  =4  or  5).  As  used  in  lUNiyriON  M’  it  must  have  unit  integral 
;uid  unit  full  width  at  half  maximum. 


Tor  example,  if  one  desires  a time  history  with  the  shape  of  the 
traj^ezoid  in  figure  4,  which  rises  in  0.25  [)ulse  length  and  stays  flat  to 
0.5  ])ulse  length,  one  would,  in  fUNCIlON'  IT,  choose 

t =0.25  t, 

1 0 

and 

t =0.5  t,. 

2 j 


Ihen  t,  and  f,  are  chosen  by 
re(|uiring  unit  full  width  at 
half  maximum. 


(t. 


t,  - t.) 


1 


;uid  unit  integral. 


2 


1 . 


These  imply 


(49) 

(50) 
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(51) 

(52) 

lliese  numbers  specify  the  trapezoidal  time  history  of  Figure  4 as  would  be 
used  in  FUNCTION  FT. 
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SECTION  6 
OUTPUT  FORMAT 


Before  the  time  cycle  printout,  an  initial  print  shows  the  fixed 
problem  parameters.  The  problem  number,  title,  and  first  input  cards  are 
jirinted  followed  by  some  of  the  XB(I),  and  SL(I),  and  all  of  KlO(l).  Tlien 
the  emission  velocities,  sjiectrum,  and  weights  are  printed,  followed  by  the 
plasma  period  (.Hquation  (42)  or  (43)),  the  parameter  a of  bciuation  (29), 

A of  liquation  (31)  (called  AL  in  the  code),  and  a ((uantity  DL  which  is 
one-half  the  particle  length  (finite  length  particles  are  used  for  the 
computation  of  charge  and  current  densities).  Finally  the  dimensional  units 
of  dynamical  quantities  arc  jirinted  in  mixed  Caussiiin  esu  and  in  MKS  units. 


Ilie  time  cycle  output  begins  with  the  cycle  number  and  the  following 


TIMI: 

NP 

,\KILL 

NPXMAX 

QTOr 

QR 


dimensionless  time  elajised 

total  number  of  live  particles 

total  number  of  particles  that  have  been  killed 
since  time  0. 

total  numT)er  of  particles  that  have  passed  XMAX 
and  been  killed. 

total  dimensionless  charge  per  unit  area  outside 
the  surface. 

total  dimensionless  charge  per  unit  area  that 
has  passed  XMAX. 


Real  time  elapsed  is  obtained  by  multiplying  "TlMi:"  by  the  unit  of  time  in 
seconds  printed  out  at  the  end  of  the  initial  j)rint.  Real  charge  per  unit 
area  is  obtained  by  multiplying  QTOT  by  the  unit  of  charge/unit  area  printed 
above,  and  similarly  for  QR . 
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Then  the  electric  field  and  potential  at  the  surface  are  printed, 
followed  by  the  field  and  potential  at  every  tenth  length  unit  from  the 
surface,  out  to  x=  10. 0.  If  fields  and/or  jiotentials  are  desired  farther 
out  tluui  this,  the  print  variable  NPR.^  may  be  changed  in  the  code.  The 
physical  field  and  potential  are  obtained  b>’  multiplying  the  dimensionless 
out]Hit  numbers  by  the  units  [)rinted  at  the  end  of  the  initial  jirintout. 

Ne.xt , the  numbers  of  particles  in  the  first  100  cells  are  printed, 
together  with  the  total  charge  (i^er  unit  area)  in  each  of  the  first  60  cells 
llien  follow  the  charge  density,  current  density  and  dipole  moment  every  0.1 
length  unit,  and  finally  the  sjiatial  integral  of  the  current  density, 

'JO 

I j(x,t)d.x  . 

0 


ilie  field  is  computed  every  cycle  at  the  edges  between  bins.  The 
field  printed  out  every  0.1  length  units  is  obtained  by  simple  linear  inter- 
polation from  tlie  fields  at  cell  edges. 


AVERAGING  FOR  CHARGE  AND  CURRENT  DENSITIES 


The  charge  density  and  current  density  are  comjiuted  every  0.1 
length  unit  by  the  following  averaging  [irocess. 


The  array  K10(  I ) shows  which  cell  contains  the  jioint  of  interest 
say  0.7  leiigth  units  so  that  I =7.  lor  ex;uni)le  in  the  space  grid  of 
Figure  .T,  0.7  length  units  is  in  cell  K10{7)  =6.  Ilic  two  adjacent  cells 
(5  and  7)  are  also  considered,  ilie  fractional  distance  y that  the  point 
x = 0.7  occui)ics  within  cell  6,  is  computed. 


x-SL(6)  _ 0.7-SI46) 

XB(6)  XB(6) 


(5.T) 
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where 


0£y<l.  (54) 

Ilien  the  current  densit)'  at  x,  j(x),  is  a weighted  average  of  the  current 
densities  C(5),  C((i)  , C(7)  in  cells  5,  6,  and  7,  respect  ivcl,v.  The  weights 
''s’  ''(■)’  given  to  cells  5,  0,  mul  7 depend  on  y in  such  a way  that  they 
have  the  following  values  when  x is  at  an  edge  (y = 0 or  1)  or  at  the 
center  of  cell  fa  (y  = 1/2)  : 


y 

0 

.5 

1 

"".s  = 

.5 

2 

0 

"fa  = 

. 5 

.0 

.5 

II 

0 

. 2 

.5 

Ihc  weights  are  given  by  the  formula 


"5 

1 M 
= 2 

(55a) 

"7 

1 PI 
= 2 >' 

(55b) 

"fa 

= 1 -w^-w. 

(55c) 

where 

>’l 

= 1.. 52 193 

(55d) 

liquations  (55)  inteqiolate  between  the  weights  given  in  the  above  table. 
When  the  point  of  interest,  x , falls  at  the  center  of  a cell  (>'=y 
averaging  gives  three  times  the  weight  to  this  cell  as  to  either  adjacent 
cell.  When  x falls  at  the  boundary  between  two  cells  (y  = 0,  or  1)  those 
two  cells  receive  ec)ual  weight. 


SECTION  7 
SAMPLE  RUN 


Wc  present  the  results  of  a samjjle  run  consisting  of  two  problems. 
Both  are  for  a 5 keV  blackboily  X-ray  spectrum  incident  on  aluminum.  The 
first  problem  has  a time  history  which  is  a symmetic  triangle  (1T=2)  with 

_ g 

full  width  at  half  maximum  of  10  ns  (FIVIIM  = I^IO  sec),  and  a fluence  of 
-4  2 

2x10  cal/cm  . The  second  problem  has  a time  history  which  is  a continu- 
ously linearly  rising  X-ray  flux  ( IT  = 1,  I’OW  = 1.0),  and  hence  there  is  no 
ITv'HM  or  fluence.  ITie  jiarameter  IT.U  in  the  second  jiroblem  is  chosen  from 

liquation  ( .X.>)  so  tliat  the  flux  is  the  same  as  the  rising  portion  of  the 

12  ’’ 

pulse  in  the  first  problem.  In  this  case,  F'LIJ  = 2 . 0x10  in  units  of  cal/cm“/ 

■>  12 
sec“.  In  both  jiroblems,  the  yield  is  Vni  = 2..S7xl0  elec/cal. 

llie  input  cards  for  this  run  are  shown  in  Figure  5.  Figure  6 shows 
the  initial  print  output  for  problem  1.0,  and  Figure  7 shows  a s;unple  of 
the  time  cycle  output.  Figure  8 shows  the  initial  ]>rint  for  problem  2.0, 
and  Figure  9 is  a sample  of  its  time  cycle  output. 

Referring  to  Figure  5,  we  have  chosen  a space  grid  out  to  x=50.0 
length  units  (essentially  50  Debye  lengths),  and  have  used  .550  cells  with  a 
spacing  ratio  7 = 1.004.  Tlie  time  step  is  0.025  time  units  (essentially 
0.025  plasma  periods)  and  we  print  every  10  cycles.  The  first  problem  will 
run  for  500  cycles,  the  second  for  400  cycles.  There  are  .52  input  cards 
defining  the  spectrum.  52  juirticles  will  be  emitted  each  time  step. 

In  Figure  6 the  initial  print  shows  the  first  four  data  cards, 
then  various  bin  sizes  XB  ;uid  distances  SI,.  I'he  full  array  KIO  is 
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SAMPLE  PPOBLEM  ONE.  5 LEV  BEl  ON  ALUMIHMM 
SOa  S50  0 0 32  2 10  1 

0.025  50.0  1.004  0.0  1.0 


2.5.rE12 

2.0E-4 

0. 15 

5.4476-06 

0 . 45 

1 .634E-05 

1 .0 

3.631E-05 

1.5 

9.7686-06 

2.0 

1 . 173E-05 

2.5 

8 . 7 74E-06 

3.0 

9. 710E-06 

4.0 

8.815E-06 

5 . 0 

6.562E-06 

b . 0 

5.35OE-0b 

r.o 

4.348E-06 

8 . 0 

3.535E-06 

10.0 

2.518E-06 

12.0 

1 .625E-06 

1 4 . 0 

1 . 048E-06 

16.0 

6 . 7bb'E— 07 

18.0 

4.390E-07 

20 . 0 

2 . 850E-0? 

22.0 

1 .571E-07 

24.0 

8 . 67yE-Li8 

26.0 

8.286E-08 

28.0 

4.252E-08 

30.0 

2.355E-08 

33.0 

1 .2356-08 

36.0 

6.527E-09 

39.0 

3.019E-09 

42.0 

2.728E-09 

46.0 

1 . 132E-09 

50.0 

4.714E-10 

55.0 

9.353E-11 

60.0 

3.502E-1 1 

ro.o 

5.813E-12 

SAMPLE  PROEiLEtl  TIJO.  LI  NEAP  PISE. 

400  350  1 0 32  1 10  0 

0.025  50.0  1.004  1.0  2.00 

2.57E12  2.E12  0. 


Figure  5.  Data  cards  for  sample  run. 
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Figure  6.  Initial  printout  for  problem  1 
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Figure  6 (Continued).  Initial  printout  for  problem  1 
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Figure  7.  Time  cycle  print  for  problem  1.0. 
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Figure  8.  Initial  printout  for  problem  2.0. 
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Figure  9.  Time  cycle  print  for  problem  2.0. 
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printed  next,  showing  for  cxainiUc  that  x = 2.4  occurs  in  cell  35  (between 
SL(35)  ajid  SL(3fa)).  Tlie  injiut  spectrum  is  printed  first  as  read  in  and 
then  after  conv'ersion  to  a noniial  energ)'  spectrum.  Hie  spectrum  is  then 
converted  to  a nonnal  velocity  spectrum  ;uid  printed  again.  The  average 
normal  velocity  is  corninited  ;ind  is  used  as  the  luiit  of  velocity.  The 
velocities  and  spectrum  are  tlien  normalized  ;uul  printed  again.  Iin;ill>’, 
equal  weights  are  cho.sen  for  the  emission  particles,  and  their  velocities 
are  printed. 

llie  plasma  period,  computed  hy  l.quat  ion  (8)  or  (IJ).  t tu'  paiMiiiet(.-r  i 
of  liquation  (29),  ;uid  the  quantity  A of  Tqu.ation  (31),  (c.illed  Al,  in  the 
code)  are  then  printed  togetlier  with  IH..  I'l.  is  one-half  t)f  tie  particle 
lengtli.  Particles  are  given  a length  of  0.  1()  length  units  (e.sentially  U.Ki 
Debye  lengths!  for  the  purpose  of  comi)uting  their  cli.irge  and  current  densities 
more  smoothly. 

Finally  the  physical  units  are  printed  in  (laussian  cgs  and  in  MRS 
units.  Quantities  printed  in  the  time  cycle  out|)ut  must  be  multijilied  by 
the  approjniate  dimensional  unit  to  obtain  the  physical  result.  lor  example, 
the  unit  vif  time  is  4.015  ns,  the  unit  of  electric  field  is  37.35  kV/m 
(or  1.246  esu) , and  the  unit  of  length  is  10.59  cm. 

As  a sample  of  the  time  cycle  outjiut , we  show  cycle  130  of  [iroblem 

1.0  in  Figure  7.  Tliis  is  at  time  3.25x4.015  ns,  or  13.05  ns.  There  are 

3364  live  particles.  A total  of  796  have  been  killed  since  the  start  of  the 

problem,  but  none  has  passed  XM/\X.  The  total  charge  is  1.507  charge  units, 

-7  2 

or  4.984x10  coulomb/m  , obtained  by  multiplying  by  the  unit  of  charge/unit 
area.  The  surface  field  is  1.507x37.35  kV/m= 56.29  kV/m. 

The  printout  for  the  fic.d  every  0.1  length  unit  shows  that  at, 
say,  0.4  length  units  from  the  surface  (4.24  cm)  the  field  strength  has 
dropped  to  0.8753x37.35  = 32.69  kV/m. 
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Hie  iiotcntial  is  measured  from  zero  at  the  surface.  M 10.59  cm 
(I  lenj^th  imit),  for  e.^ample,  it  has  risen  to  0. 853.5X.5.955  kV  = 3.375  kV  as 
shown  in  the  printout. 

llic  next  group  of  numbers  shows  how  many  particles  are  in  each 
cell.  There  are  285  in  the  first  cell,  221  in  the  second,  etc.  The  charge 
(per  unit  area)  in  the  first  several  bins  is  printed  next. 

The  charge  density  and  current  density  every  0.1  length  unit  are 
printed  next  in  dimensionless  form.  The  first  number  printed  is  at  0.1 
length  unit  (1.059  cm)  from  the  surface,  not  at  the  surface.  The  surface 
charge  and  current  densities  are  not  computed  by  the  code.  Thus  at  130  cycles, 
the  charge  density  at  1.059  cm  is  2.015x5.122x10  ^ = 6.29x10  ^ Coul/m'^  , 
and  the  current  density  is  0.1369x82.35  = 11.27  Am])s/m“. 

The  electron  number  density  can  also  be  obtained  from  the  dimen- 
sionless charge  density  outj)ut  by  multipl\ing  by  the  unit  of  number  density 
7 3 

(1.949x10  elec/cm"),  showing,  for  example,  that  at  0.6  length  units  (6.35  cm) 

*7  7 _ X 

the  number  density  is  0.6729x1.949x10^  = 1.31x10  cm  '. 

llie  dipole  moment  (per  unit  surface  area)  contributed  by  all 
electrons  out  to  x, 

X 

I’  = ^ xp(x)  dx 
0 

is  printed  next.  Mere  the  first  number  is  at  the  surface,  x = 0,  and  so 
is  0.  At,  say,  0.7  length  units  (7.41  cm)  the  dipole  moment  is  seen  to 
be  0.2212x3.502x10  ^ = 7.75x10  ^ Coul/m.  Out  to  9.9  length  units  it  is 
1.511x3.502x10*^  = 5.29X10*^  Coul/m. 
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linall)-,  the  integral  of  the  current  densits-, 


t:i  = ^ j (x,t)  dx 
0 

is  printed,  having  the  value  0.8625x8.722  = 7.52  Amps/m. 

The  initial  print  for  problem  2.0  is  shown  in  Figure  8.  The  Xh , 

SI,,  and  KIO  arrays  are  the  same,  and  the  spectrum  is  the  same.  The  plasma 
period  in  this  case  is  computed  by  liquation  [45)  w i tli  I’OW  = 1.0,  and  h.is  a 
slightly  different  value  than  in  problem  1.0.  Similarly,  the  units  of 
length,  number  densit)’,  etc.  are  not  exactly  the  same. 

We  show  the  output  for  cycle  50  in  Figure  9.  It  is  at  time 

1.25  units  = (>.804  ns.  llie  fastest  emitted  particle  had  velocity  5.827  units 
and  so  will  not  have  travelled  farther  than  4.78  length  units  by  this  time. 
Hence,  the  field  beyond  4.7  length  luiits  is  still  0,  ;uid  the  printout  shows 
no  particles  beyond  the  64-th  cell.  Otherwise  the  interpretation  of  quantities 
is  the  same  as  for  the  previous  problem. 

Problem  1.0  ran  for  500  cycles  ;uid  moved  a maximum  of  5542  jiarticles 
(at  cycle  190).  Problem  2.0  ran  for  400  cycles  :ind  moved  a maximum  of  2500 
particles  fat  cycle  110).  Total  Cl’  time  for  both  problems  on  a CDC  7600 
was  41.4  sec. 


44 


DISTRIBUTION  LIST 


DKF’ARTMKNT  OV  DKKfiNSi: 

Director 

Detense  Advanced  Researclj  I’roj.  Agency 
ATTN ; N'MR 

Director 

Detense  (i<itmnun icat  Ions  Agency 
ATT..;  NMR 

Delense  lX)c\imentat  ion  Center 
(^meron  Station 
12  cy  ATTN:  TC 

Diret  tor 

Defense  Intelligence  Agency 
ATTN:  DB-^^C 


Director 

De tense  Nuclear  Aj’.ency 
2 cy  ATTN:  R.\KV 

ATTN:  STSl,  Archives 

} cy  ATTN:  STTL,  Tech.  Library 

ATTN:  DDST 


Cot;  inlander 
Kield  Command 
Defense  Nuclear  Agency 
ATTN:  KCLMC 

ATTN-. 


Director 

Interservice  Nuclear  Weapons  School 
ATTN:  Document  Control 


Director 

Joint  Strat.  Target  Planning  Staft,  JCS 
ATTN:  Jl.TW-2 


Chief 

Livermore  Division,  Field  Comnuind,  DNA 
Liiwrence  Llvernw>re  Uihoratory 
ATTN:  FCPRL 

National  Communications  System 
Office  of  the  Manager 
ATTN:  NCS-TS 


Director 

National  Security  Agency 
ATTN:  K-42^ 

OJCS/J- j 

ATTN;  J-3,  RDTA  Br.,  WWMCCS  Plans  lUv. 
OJCS/J-5 

ATTN;  J”5,  Plans  & Policy  Nuc.  Dlv. 

Dir.  of  Defense  Research  & Knglneering 
Department  of  Defense 
ATTN:  S&SS  (OS) 

DKPAKTMKNT  OF  THK  ARMY 

Director 

HMD  Advanced  Tech.  Center 
Huntsville  Office 

ATTN;  RDMH-0 


DKPARTMFNT  OF  THK  ARMY  (Cotuinued) 

Comitiaiuler 

BMD  System  Command 

ATTN:  BDMSC-TKN 

Dep.  Chief  of  Staff  for  Hesearcli  Dev.  & Acq, 

De[>artment  of  the  Army 
ATTN:  DAMA-CSM-N 

Commander 

Harry  Dl.imond  Laboratories 

ATTN:  DKXDO-RCC,  John  A.  Rosado 

ATTN:  DRXDO-TI,  Tech.  Lib. 

ATTN:  DRXDO-RCC,  Raine  Gilbert 

Commander 

Pleat Inny  Arsenal 
ATTN:  SARPA 

ATTN:  SMIT'A 

Commander 

Redstone  Scientific  Information  Center 

US  Army  Missile  Command 

ATTN:  Chief,  Documents 


Chief 

US  Army  Communications  Sys.  Agency 
ATTN : SCCM' AD- SV , Library 


Commander 

US  Army  Klectronics  Command 
ATTN:  DRSKL 

Commander 

US  Army  Foreign  Science  & Tech.  (!enter 
ATTN;  DRXST-ISI 

DIIPAKTMKNT  OF  THK  N>VY 

Chief  of  Naval  Operations 
Navy  Department 

ATTN:  Code  bfViC) 

Chief  of  Nava]  Research 
Navy  Department 

ATTN:  Henry  MuIIaney,  Code  ^27 

Director 

Naval  Research  Laboratory 
ATTN;  Code  7701 
ATTN:  (U>de  7 756,  .1,  Davis 


Commander 

Naval  Surface  Weapons  Center 

ATTN!  Code  WA501 , Navy  Nuc.  Prgms.  Off. 

Director 

Strategic  Systems  Project  Office 
Navy  Department 
ATTN;  NSP 


45 


DH’AHfMKNr  OF  THF  AIK  FOKCK 


DFPARTMFNi;  OF  DKFFNSF  CONTRACTORS  (Cont  imit'tj) 


AK  C.fophysics  t-ah»>ratory  ♦ AFSC 
ATTN:  t.lKirli's  Pikt* 

AF  M.»tfrials  Laboratory.  AFSC 
ATTN:  Library 


AF  We.ipi>ns  1 

.abor.itory 

2 cy  ATTN: 

DYC 

AI TN: 

NT 

ATTN : 

SUL 

2 cy  ATTN: 

NTS 

Hq.  USAF/RD 

ATTN: 

RDQSM 

Conuiuindcr 

Homo  Air  Developmont  Contor,  AFSC 
ATTN:  Kdwarcl  A.  Biirko 

S/\MS0/IJY 

ATTN:  DYS 

SiVMSO/MN 

ATTN:  MNNH 

ATTN:  MNNC 

S/\MSO/SK 

ATTN:  SKK 

SAMSO/XR 

ATTN : XRS 

tommandor  In  Cliiof 

Slratfj^lc  Air  Comiruuul 

ATTN:  NRI-STINFO.  Library 

ATTN:  XFFS 

ENLRGY  RKSKARCH  & DFVKLOFMFNT  ADMINISTRATION 


Inlversity  of  California 
Lawrence*  Livermore  Laboratory 
ATTN:  l.-yn 

ATTN:  Into.  Dept..  L-'i 

Ixis  Alamos  .Scientific  Laboratory 

ATTN:  Doc.  Contrt>l  for  Report.s  Lib. 

Sand  la  Liboratories 
Livermore  Laboratory 

ATTN:  lior.  Con.  for  Theodore  A.  Dell  in 

Sandia  Laboratories 

ATTN:  iJoc . Con.  for  3141  Sandia  Rpt . Coll. 

0Tllt.k  COVFKNMKNT  A(;KNC1ES 
NASA 

Levis  Research  Center 

ATTN:  N.  .1.  Stevens 

ATTN:  Carolyn  Purvis 

ATTN:  Library 

UKPARTMLNT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corporation 
ATTN:  l.ihrary 

ATTN:  V.  Josephson 

ATTN:  Frank  Hal 

ATTN:  Julian  Reinhoimer 


Avco  Research  4 Systt*ms  Croup 

ATTN:  Research  Mb..  AHiO.  Km.  7201 

’nie  Boeing  Comp.any 

ATTN:  Preston  (iert*n 

University  of  California  at  San  Diego 
ATTN:  Sherman  De  Forest 

C<>mputer  Sciences  Corporation 
ATTN:  Alvin  T.  Schiff 

Dr.  Eugene  P.  DePloml) 

ATTN:  Eugene  P.  DePlomb 

Dikewood  Industries,  Inc. 

ATTN:  Tech.  Lib. 

ATTN:  K.  Lee 

E(;6.C.  Inc. 

Albuquerque  Division 

ATTN:  Technical  I.ibrary 

Ford  Aerospace  & Communications  Corp. 

ATTN:  Library 

ATTN:  Donald  R.  McMorrow,  MS  C30 

Cen e ra 1 E 1 ec trie  Company 

Space  Division 

Valley  Forge  Space  Center 

ATTN:  Joseph  C.  Peden,  VFSC.  Rm.  4230M 

(ieneral  Electric  Company 
TEMPO-Center  for  Advanced  Studies 
ATTN:  William  McNaraina 

ATTN:  DAS 1 AC 

Hughes  Aircraft  Company 
Centincla  4 Teale 

ATTN:  Tech.  Lib. 

Hughes  Aircraft  Company,  El  Segundo  Site 
ATTN:  Edward  C.  Smith,  MS  A620 

ATTN:  William  W.  Scott,  MS  AlORO 

Institute  for  Defense  Analyses 
ATTN;  IDA  Librarian 

IRT  Corporation 

ATTN:  Dennis  Swift 

ATTN:  Library 

Jaycor 

ATTN:  Eric  P.  Wenaas 

ATTN;  Library 

Jaycor 

ATTN:  Robert  Sullivan 

Johns  liopkin.s  University 
Applied  Physics  Laboratory 

ATTN:  Peter  E.  Partridge 

Kaman  Sciences  Corporation 

ATTN:  Jerry  I.  Lubell 

ATTN:  W.  Foster  Rich 

ATTN:  Library 


46 


DF.PARTMFNT  0F_  UFFKNSE  .(’ON  rRACTORS  j^l'^l'lNSJ^  CONTWCTOJ^S  (Cont  Inued) 

Ix>ckht't'ti  Missiles  & Space  Company,  Inc.  Science  Applications,  Incorporated 

ATTN:  Dept.  85-85  ATTN:  Will  lam  L.  Chadsey 


McDonnell  Douglas  Corporation 
ATTN:  Stanley  Schneider 


Simulation  Physics,  Inc. 

ATTN:  Roger  Little 


Mission  Kesea 
ATTN 
ATTN 
5 cy  ATTN 
ATTN 


rch  Corporation 
Roger  Stettner 
Conrad  L.  Longmire 
Library 
N.  J.  Carron 


Mission  Research  Corporal ion-San  Diego 
ATTN:  V.  A.  J.  Van  Lint 

ATTN:  Library 

K & D Associates 

ATTN:  Leonard  Schlessinger 

ATTN:  Library 

Rockwell  International  Corporation 
ATTN:  N.  J.  Kudie,  FA53 

ATTN:  Library 


Stanford  Research  Institute 
ATTN:  Library 

Systems,  Science  & Software,  Inc. 
ATTN:  Andrew  R.  Wilson 

ATTN:  Library 


TRW  Systems  Croup 

2 ry  ATTN:  Robert  M.  Webb,  Rl-2410 

ATTN:  Tech.  Info.  Center/S-1 930 


47 


